Non-Abelian geometric phases in ground state Josephson devices 
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We present a superconducting circuit in which non-Abelian geometric transformations can be 
realized using an adiabatic parameter cycle. In contrast to previous proposals, we employ quantum 
evolution in the ground state. We propose an experiment in which the transition from non-Abelian 
to Abelian cycles can be observed by measuring the pumped charge as a function of the period of 
the cycle. Alternatively, the non-Abelian phase can be detected using a single-electron transistor 
working as a charge sensor. 

PACS numbers: 03.65. Vf, 85.25.Dq 



Introduction — Accurate control and measurement of 
few-level quantum systems has recently attracted great 
experimental and theoretical interest with possible appli- 
cations in quantum information processing (QIP). Geo- 
metric phases^ arising from adiabatic and cyclic quantum 
evolution can provide robustness against, e.g., timing er- 
rors. Recently, it was shown that such evolution in a 
non-degenerate ground state is immune to decoherence 
from a low-temperature environment 2 suggesting that it 
may provide an important tool for controlling quantum 
systems. 

In the non-degenerate case, the accumulated geomet- 
ric phase, the Berry phased, is a shift of the complex 
phase of the eigenstate, and hence cannot be used as 
such for QIP. Non-Abelian phases^ correspond to uni- 
tary matrices operating in a degenerate subspace of the 
system Hamiltonian, thus providing means for universal 
QIP 6 . Although geometric phases capable of entangling 
two quantum bits, qubits, have been observed in liquid- 
state nuclear magnetic resonance experiments^, this kind 
of geometric quantum computing (GQC) is yet to be 
demonstrated. In fact, the geometric phases using nu- 
clear magnetic resonance^, and in more recent experi- 
ments^ demonstrating non-adiabatic Aharonov-Anandan 
phases^ in fullerene spin qubits, accumulate in a rotat- 
ing frame, and hence there is no true degeneracy in the 
system. 

The initial proposals for the experimental realization 
of GQ C 10 i n rely on a fully degenerate subspace to build 
the logical operators and it has been extended to many 



quantum systems* 2 -— including Josephson junction de- 
vice s] 15 ' 16 . In similar systems, a way to observe the 
non-Abelian evolution by measuring the charge pumped 
through the device has been recently proposed^. 

However, all the schemes assume typically a so-called 
tripod Hamiltonian which has degeneracy only in its ex- 
cited states [see Fig. QJa)] rendering the system prone 
to decoherence even in the low-temperature limit. This 
is potentially a serious limitation in the condensed mat- 
ter systems in which the coupling between system and 
environment is strong and unavoidable. 

In this paper, we present an experimentally realizable 
Josephson device and show that it can be used to observe 
adiabatic non-Abelian geometric phases. In contrast to 
the above pioneering works, we employ a conceptually 
different Hamiltonian allowing us to work on the ground 
state manifold of the system. This proposal provides a 
clear extension to the theoretical proposals^— and ex- 
periments 2 ^' 2 ^ on the Berry phase in superconducting cir- 
cuits. 

Non-Abelian adiabatic evolution — We denote the pa- 
rameters of the system Hamiltonian in a general cyclic 
loop by a vector A. The instantaneous eigenstates of 
the Hamiltonian H[X(t)] along this loop for all t € [0, T] 
are denoted by {\ip a (t))}, where T is the period of the 
cycle. Generally, any temporal evolution of the system 
state can be represented using the time evolution opera- 
tor, U(T), such that |*(T)) = [/(T)|*(0)), where 
is the state of the system at time t. The charge Q trans- 
ferred through a superconducting system in one parame- 
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ter cycle can be obtained by integration of the current op- 
erator I = -^d v H(t) as Q = J T dt (*(*)|I|#(t)), where 
and tp is the superconducting phase difference across the 
syste m 17 ! 22 . Using the Schrodinger equation and the def- 
inition of the time evolution operator, this can be written 
in the form 



Q = -2ia{*(0)\U\T)ld v U(T)}\*(0)). (1) 

However, if the Hamiltonian parameters are changed 
adiabatically along the cycle, the evolution can be re- 
stricted to the initial eigenspace. In an n-fold degener- 
ate eigenspace, the state of the system after a param- 
eter cycle is |*(T)) = C/ ad (T)|*(0)) + <D(1/T), where 
I^W) = J2^=i I If the instantaneous eigen- 
vectors are defined globally and continuously, the opera- 
tor £/ad(t) is represented in this basis as 



(2) 



where E(t) is the energy of the degenerate eigenspace, T 
is the time ordering operator, and the connection T(t) is 
given by [r(t)] a/9 = (ip a (t)\ipp(t)). The first exponential 
function in Eq. (J2|) yields the accumulated dynamic phase 
shift, Ud yn {t), and the second one provides the geometric 
transformation, U geo (t), which is non-Abelian in general. 

In the adiabatic limit, U(T) can be replaced with 
U ad (T), and substituting Eq. © into Eq. ((TJ) yields the 
relation between the different transformations and the 
transferred charge 

Q = -2ie((*(0)|E/+ eo (T) [d v U gco {T)] |tt(0)> 

+ (*(0)|C/t yn (T) [d v U dyn (T)} |*(0)>), (3) 

where the first term is the geometric pumped charge 
and the second the dynamic charge due to the usual su- 
percurrent. In the case of a nondegenerate eigenspace, 
n = 1, this reduces to the well known relation Q — 
2ed v (&B — < dd), where the accumulated Berry phase, 0b, 
is related to U geo by U gco = e l&B and the dynamic phase, 
0d, to Udyn by Udyn — e _l6ld ^. See Brosco et al*L for an 
alternative way to obtain the pumped charge. Although 
the Berry phase induces just a phase shift to the state 
vector, it does not commute in general with the current 
operator I which originates from a higher-dimensional 
system. 

Model circuit — The Cooper pair pump shown in 
Fig. [ITc) is considered here as the physical realization 
for observing non-Abelian geometric phases. It consists 
of three SQUIDs in series with two superconducting is- 
lands between them. The SQUIDs are operated as tun- 
able Josephson junctions which can be closed (Josephson 
energy Ei is zero) and opened (Ei ^ 0) by controlling the 
magnetic flux through them. The phase difference of the 
order parameter across the whole device, tp — c£>l — (f>R, is 



kept constant by the magnetic flux $ through the outer- 
most loop. The Hamiltonian has five external parameters 
which are controlled during a pumping cycle, i.e., three 
magnetic fluxes and two gate voltages. 




FIG. 1: (color online) (a) Eigenenergy diagram of the so- 
called tripod Hamiltonian consisting of two bright states, B, 
and two degenerate dark states, D. (b) Eigenenergy diagram 
of the circuit Hamiltonian along the considered cycle. The 
energy difference between the ground, g, and excited, e, state 
is denoted by A and the ground state degeneracy splitting by 
8. In the ideal case, 8 = 0. (c) Schematic diagram of the 
Non-Abelian superconducting pump. Green and blue denote 
the superconducting islands and red the Josephson junctions. 

The charging energy part of the Hamiltonian, i? c h, is 
given by 



H ch =E Cl (ni - n gl ) 2 + E C ' 2 (n 2 - n g2 ) 2 
+E m (hi - n g i)(n 2 - n g2 ), 



(4) 



where hi is the operator for the excess number of Cooper 
pairs on the ith island and n gi is the corresponding gate 
charge given by n g i — C g iV g i / (2e). The charging energies 
are E Cl = 2e 2 C J2 JC 2 , E C2 = 2e 2 C7 Ei /C7 2 , and E m = 
Ae 2 C m /C 2 . Here, C E . is the total capacitance of the ith 
island and C 2 = C El C7 E2 - C 2 m . 

The Josephson part of the Hamiltonian, Hj, reads 



Hj= ^2 {j c ss\ni + l,n 2 ){n 1 ,n 2 \ 
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+Jcff,mK + l,n 2 - I)(ni,n 2 | 
+JeflF,2|ni,na + l)(ni,n 2 | +h.c.J, (5) 

where \ni,n 2 ) denotes the state with rii excess Cooper 
pairs on the ith island, J cS ,i = -^i($i)e 4v(<E, ^ 2 /2, 
Joff, 2 = -S 2 ($ 2 )e-^(*)/ 2 /2, and J eff , m = -£ m ($ m )/2. 
Here, E\ , E 2 , and E m are the tunable Josephson energies. 
The full Hamiltonian is given by H = i? c h(F^i, V g2 ) + 
ff. 7 ($i,$ 2 ,$ m ,$). 

Non-Abelian cycle — If all the SQUIDs are closed, the 
conventional stability diagram with a hexagonal struc- 
ture is recovered 2 ^, see Fig. [2] In the vicinity of the triple 
degeneracy point of states |1,0), |0, 1), and |f,I), the 
adiabatic evolution is approximately restricted to these 
three states. The parameter cycle is composed of three 
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symmetric paths in each of which a SQUID is opened, the 
gate voltages are shifted along a ground state degeneracy, 
and finally the SQUID is closed. 

Along each path, the effective three-level Hamilto- 
nian has a 2 x 2 block and can be written as -ff e ff = 



a 



B(t) + ek(t)\k)(k\, where <?ij = {a. 
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a vector composed of the Pauli matrices for the states 
i,j (for example, of . = + \j){i\), B(t) is an effec- 

tive magnetic field, e& is the eigenvalue of the third charge 
state, and{|i),|j),|fc)} = {|l,0),|0,l),|l,l)}. 

The condition of the ground state double degeneracy 
is satisfied by tuning the smaller eigenvalue of the 2x2 
block of -ff e ff to be equal to e/-(i) along the evolution. In 
the three-level approximation this implies that the de- 
generate gate voltage paths are hyperbolas in the gate 
voltage plane with one SQUID kept open. Along the 
opening and closing of the SQUIDs, we choose to change 
voltages linearly with the SQUID energies. In this way, 
a nontrivial loop encircling the triple degeneracy point 
can be traversed along a path with a doubly degenerate 
ground state. 

Using the eigenstates along the three paths, we can 
construct a continuous global basis (defined in the whole 
parameter space) and calculate the connection [r(t)] a ^. 
If the SQUIDs can be perfectly closed, the supercurrent 
due to the dynamic phase in Eq. ([3]) vanishes since the 
energies of the eigenstates do not depend on (p. In this 
case, the transferred charge has only a geometric contri- 
bution which can be calculated analytically from the tp 
dependence of the U gco (T) operator. For a cycle starting 
from the degeneracy line between the states 1 1 , 0} and 
0, 1), this yields for the geometric transformation 



U geo (T) = 



e l f 

1 



(6) 



represented in the basis {|1, 0), |0, 1)}. This result was 
confirmed by solving numerically the Schrodinger equa- 
tion using 25 charge states indicating that our analysis 
does not rely on the three-state approximation. The ob- 
tained transformation is topological in the sense that it 
does not depend on the exact values to which the SQUIDs 
are opened as long as the evolution is kept degenerate 
along the cycle. From Eq. ([3]), we obtain for the geomet- 
rically pumped charge 



Q = 2e[c;(0),c5(0)] 
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2e|c 2 (0)| 2 . (7) 



Thus, the pumped geometric charge is independent of the 
phase across the device and depends only on the initial 
state. 

Observation scheme for non-Abelian phases — Here we 
discuss two methods to observe the non-Abelian trans- 
formations. Firstly, a single-electron transistor (SET) 
can be coupled asymmetrically to the superconducting 
islands and used as a charge sensor. The additional ca- 
pacitance due to the SET changes slightly the charging 



FIG. 2: (color online) Parameter cycle inducing non-Abelian 
transformations with the stability diagram as the background. 
The charge state \i,j) which minimizes the charging energy in 
each hexagonal area is denoted by The z axis represents 
Ei with blue in the front for E\, purple on the right for E m , 
and red on the left for Ei. The SQUID energies can be closed 
down to -Ef m , which is zero in the ideal cycle. The inset shows 
the parameter cycle projected into the gate voltage plane. 



energies but does not affect the operation principle of the 
circuit. Initializing the system and performing the pa- 
rameter cycle adiabatically swaps the charge states of the 
islands regardless of the phase across the device, which 
can be detected with the charge sensor. Observation of 
this charge transfer proves the non-Abelian character of 
the evolution since in the Abelian case, initial popula- 
tions are conserved in cyclic adiabatic evolutions. 

Another way to observe the non-Abelian features is 
to measure the pumped charge through the system us- 
ing a detector junction 22 . Since in the experiments the 
SQUIDs cannot be perfectly close d 15 ' 24 , we consider here 
a case in which the Josephson energies can be tuned down 
to Ef 11 " 1 of their maximum value E™ ax . 

In the case of non-ideal SQUIDs, two additional effects 
have to be considered. Firstly, the supercurrent contribu- 
tion usually dominates over the geometric contribution. 
However, it has been showr* 2 ^ that the supercurrent con- 
tribution can be efficiently measured by traversing the 
parameter cycle first forwards and then backwards. In 
the perfect adiabatic limit, the geometric component of 
the current cancels itself and the measured total current 
is twice the supercurrent. 

Secondly, the two lowest-energy eigenstates are not 
perfectly degenerate with the energy gap S ~ £™ m , see 
Fig.QJb). To obtain non-Abelian evolution, the loop has 
to be traversed fast enough such that the two lowest 
eigenstates are effectively degenerate, that is T <C h/6 
where T is the cycle period. On the other hand, the evo- 
lution should be slow enough to avoid transitions to the 
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higher states implying T ^> ft/ A, where A is the energy 
gap to the excited state. To obtain the Abelian limit, 
the energy gap S can be increased by larger Josephson 
energies and the cycle can be traversed slower such that 
no transitions occur. 

The system can be initialized to the state [1,0) by the 
following procedure. First, all the SQUIDs are closed to 
Ef 1 ™ and gate voltages tuned to have |1,0) as a nonde- 
generate ground state. After the system has relaxed to 
the ground state, the gate voltages are suddenly shifted, 
Tshift *C ft/<5, to the degeneracy line between the states 
1,0) and 0, 1). The sudden shift keeps the system in the 
state 1 1,0) and the non- Abelian cycle can be traversed 
starting from a well-known initial state. The system can 
be initialized to the state |0, 1) with a similar procedure. 

To describe the adiabaticity of the evolution, we intro- 
duce the adiabaticity parameter a defined as the popu- 
lation of the initial state after a back-and-forth cycle. In 
the perfectly non-Abelian regime, the geometric trans- 
formations induced by the forward and backward cycles 
exactly cancel each other. Thus, the total transforma- 
tion is proportional to the identity implying that a = 1. 
For the perfectly Abelian limit, no transitions occur be- 
tween the eigenstates and again a = 1 if the initial state 
is an eigenstate. Between these two regimes, no easy the- 
oretical prediction can be made since the states are only 
partially mixed during the evolution. 

In all numerical simulations, we fix the phase across 
the device (p to zero and Ec t = 0.2 meV. Figure OUJa) 
suggests that for non-Abelian cycle with period 5 ns 
< T < 10 ns the evolution is adiabatic and a is close 
to unity even if the SQUIDs cannot be perfectly closed 
with k = Ef^jEf^ = 1000. In this regime, the pumped 
charge shown in Fig. EJb) reaches the value 2e or 
depending on the initial state as predicted by Eq. (J7]). 
With k = 5000 the adiabatic evolution window is broad 
and observed as a pumped charge plateu. To obtain a 
measurable current with a reasonable averaging time (> 
1 pA)2£, the pumping cycle needs to be repeated fast 
enough. If simply a sequence of repeated pumping cycles 
is performed, the measured current reflects the average 
pumped charge e regardless of the initial state due to 
the swapping between the states |0, 1) and |1,0). On the 
contrary, the system can be initialized to the same state 
before every cycle. In this case, the pumped charge per 
cycle is 2e or depending on the initial state. Measuring 
such dependence on the initialization indicates that the 
charge states are swapped after each cycle providing a 
fingerprint of the non-Abelian evolution. 

The evolution can be made Abelian by increasing the 
cycle period and keeping all the SQUIDs constantly open 
with ^ max = Ef m = -0AE C - Figure l^c) indicates that 
the evolution is adiabatic with cycle periods longer than 
^85 ns and additional results (not shown here) confirm 
that the evolution is Abelian. Numerical simulations for 
the pumped charge, shown in Fig. EJd), yield 8.5e or 
depending on the initial state which are the two lowest 
eigenstates. The pumped charge in the Abelian limit 



depends on tp which in the simulation is fixed to zero. 
Here, the two different procedures (with and without ini- 
tialization between the pumping cycles) lead to the same 
average pumped charge pointing out that no swapping 
between the lowest eigenstates takes place. 
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FIG. 3: (color online) (a) Adiabaticity as a function of the 
cycle period in a cycle with £" lax = -0AE c and £™ in = 
Ei^/k. Red (blue) denotes a cycle starting from the initial 
state 1 1,0) (|0, 1)) for k = 1000 (dots) and k = 5000 (dashed 
line), (b) The pumped charges corresponding to (a), (c) Adia- 
baticity of the ground (red) and excited state (blue) in a cycle 
with all the SQUID energies fixed to £f ax = Bf in = ~0AE c . 
(d) The pumped charges corresponding to (c). In (b) and (d) 
the gray lines denote the geometric pumped charges in the 
perfectly adiabatic limits. The shaded areas indicate the cycle 
periods with which the evolution is adiabatic. The charging 
energy Ed used is experimentally realizable 0.2 meV— and 
the phase across the device is fixed to zero, tp = 0. 



In conclusion, we have presented a rather simple su- 
perconducting circuit with which non-Abelian geometric 
transformations can be realized in the ground state of the 
system. A parameter cycle is introduced for which the 
corresponding geometric transformation is determined 
analytically. Two observation schemes are presented for 
the non-Abelian features taking into account the most 
important experimental restrictions. 
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